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We have carried out the largest search for stable particles with fractional electric charge, based 
on an oil drop method that incorporates a horizontal electric field and upward air flow. No evidence 
for such particles was found, giving a 95% confidence level upper limit of 1.17 x 10~ 22 particles per 
nucleon on the abundance of fractional charge particles in silicone oil for 0.18e < \Q T esiduai\ < 0.82e. 
Since this is the first use of this new method we describe the advantages and limitations of the 
method. 

PACS numbers: 14.80.-j, 13.40.Em, 14.65.-q, 47.60.+i 



I. INTRODUCTION 



We have carried out the largest search for fractional 
electric charge elementary particles in bulk matter us- 
ing 70.1 mg of silicone oil. That is, we looked for stable 
particles whose charge Q deviates from Ne where N is 
an integer, including zero, and e is the magnitude of the 
charge on the electron. No evidence for such particles 
was found in this amount of silicone oil. We used our 
new version of the Millikan oil drop method contain- 
ing two innovations compared to the classical method 
that we used in Halyo et al.^. One innovation is that 
the drop charge is obtained by observing the drop mo- 
tion in a horizontal, alternating electric field compared to 
the classical use of a vertical electric field g, § f§. The 
other innovation is the use of an upward flow of air to 
reduce the vertical terminal velocity of the drop, which 
enabled us to use larger drops, about 20.6 /im in diam- 
eter compared to the 10 /im drops used in our previous 
experiments. 

We define the residual drop charge, Q r — Q — JV/e 
where Ni is the largest integer less than Q/e. We find 



the 95% confidence level upper limit on the abundance 
of fractional charge particles in silicone oil for 0.18e < 
Q r < 0.82e is 1.17 x 10~ 22 particles per nucleon. This 
experiment was a follow up on our previous search in 
silicone oil, Halyo et aZ.[|), based on 17.4 mg. In that 
search we found one drop with anomalous charge, but no 
such charge was found in the present experiment. 



In this paper we describe t he e xperimental method and 
apparatus in Sec. II. In Sec. Ill we discuss the measure- 



ment precision resulting from the various measurement 
errors and the calibration methods. The data analysis 
method, including the criteria used to accept drop charge 
measurements, is discussed in Sec. IV. Here we pay par- 



*Work supported by Department of Energy contract DE-AC03- 
76SF00515. 



ticular attention to the drop spacing criterion necessi- 
tated by interactions between adjacent drops. This is 
the primary limitation on the rate at which drops can be 
measured and we had to acquire considerable experience 
to understand this limitation. We conclude with Sec. [v|, 
giving our results, comparing our results with other frac- 
tional charge searches, and discussing the applicability 
and extension of this new Millikan oil drop technique to 
other searches. 
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FIG. 1: Basic principles of the experimental method. 



II. EXPERIMENTAL METHOD AND 
APPARATUS 

A. Experimental method 

The principle of the experimental method is simple. 
Consider a drop of radius r, density p, and charge Q 
falling in air through a horizontal electric field of strength 
E, as shown in Fig. [l[ Applying Stokes' law the horizon- 
tal terminal velocity v x is 



QE 
Q'K'qr 



(1) 



where r\ is the viscosity of air. Hence measurin g v x g ives 
Q providing r is known. As explained in Sec. Ill A , the 
drop radius is determined from the v x of integer charge 
particles. Note that the measurement of Q does not de- 
pend on the density of the drop and is also independent 
of the gravitational force on the drop. The electric field 
alternates in the +x and — x direction so that the drop 
is moved back and forth along the x axis. This cancels 
some sources of error and allows the drop motion to be 
viewed in a relatively narrow horizontal area. Previous 
uses of a static horizontal electric field in the Millikan oil 
drop method were in 1941 by Hopper and Laby S who 
measured the electron charge and by Kunkel p] in 1950 
who measured the charge on dust particles. 

If the drop were falling in still air, the vertical terminal 
velocity would be given by 



9?7 



(2) 



Hence the net downward velocity of the drop is 

9r] 



(3) 



As explained in the next section we want v z to be 
small hence we set v a i r to be close to v z , term but slightly 
smaller. 



B. General description of experiment 

Figure |H is a schematic picture of the apparatus. Drops 
averaging 20.6 fjxa in diameter are produced at a rate of 
1 Hz using a piezoelectrically actuated drop-on-demand 
microdrop ejector. The drops fall through the upward 
moving air in the measurement chamber passing through 
a horizontal, uniform, alternating electric field. In this 
figure the electric field is perpendicular to the paper. The 
electric field alternates as a square wave with a frequency 
of 2.5 Hz and has an amplitude of about 1.8 x 10 6 V/m. 

A rectangular measurement region 2.29 mm in the x 
direction by 3.05 mm in the z direction is projected by 
a lens onto the charge-coupled device (CCD) sensor of a 
monochrome, digital video camera. A light source con- 
sisting of a bank of light emitting diodes (LEDs) provides 
10 Hz stroboscopic illumination. As the motion of the 
drop carries it through the measurement region, its im- 
age appears on the surface of the CCD. Thus the camera 
collects 10 frames per second, the drop appearing as a 
dark image on a bright background. 

In addition to the v x motion there is also the v z motion, 
Eq. (^). Since the camera has a field of view of Z = 
3.05 mm in the vertical direction, the 10 Hz stroboscopic 
illumination leads to acquisition of 



10Z/v z = 30.5/u* 



(4) 



where g is the acceleration of gravity. However we use 
an upward flow of air of velocity v a i r in the — z direction. 



images of any given drop, before the drop moves below 
the viewing area of the camera. Here v z is in mm/s. 
Hence we get a larger number of images per drop, leading 
to better charge measurement precision, when v z is small. 
Of the order of Ni mages — 15 are required. 

We give an example of the importance of the upward 
airflow in obtaining this many images. Consider a typical 
drop of diameter 20.6 /mi with a density of 0.913 g/cm 3 . 
From Eq. (^|), v z .t er m = H-3 mm/s. If there were no 
upward airflow there would be an average of 2.6 images 
per drop. To obtain Ni mages = 15, v z must be about 
2.0 mm/s. Therefore from Eq. (|^), v a i r must be 11.3 — 
2.0 = 9.3 mm/s. 

Each image from the CCD camera is processed though 
a framegrabber in a conventional desktop computer, the 
signal in each pixel being recorded. An analysis program 
then finds the drop images and calculates the x and z 
coordinates of the centroid of the drop image. Using all 
the images of the drop and knowing the time spacing of 
the images, v x and v z are then calculated. 
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FIG. 3: Drop generator. 



C. Drop generator 

The drop generator, Fig. ||, is based on the general 
design principles used in piezoelectrically actuated, drop- 
on-demand, inkjet print heads. Our generator is designed 



for flexibility, allowing a variety of liquids to be used and 
providing ease of control and maintenance. The body of 
the generator is glass so as to preserve the purity of the 
liquid, but the ejection aperture at the bottom is micro- 
machined silicon 0. The diameter of the aperture sets 
approximately the diameter of the drop. Upon applica- 
tion of a short voltage pulse, usually 3 to 20 fis, across the 
surfaces of the piezoelectric disk, the central hole in the 
disk contracts in diameter, squeezing the glass tube and 
sending a pressure pulse down the liquid, ejecting a jet of 
fluid from the aperture. The forming of a discrete fluid 
drop from the high speed jet is a complex process with 
the repeatability of the process and the final diameter of 
the drop being highly dependent on the properties of the 
fluid, and on how the fluid is driven. For this reason, the 
shape and amplitude of the voltage pulse applied to the 
piezoelectric disk must be specifically tuned for stability 
and the desired drop size. In addition, it is necessary to 
experiment with single and double pulsing, varying both 
the pulse width and the separation between the pulses. 

The pressure in the drop generator is maintained 
slightly below atmospheric pressure by 10 to 30 mm Hg. 
This helps to retract the excess liquid outside the ejec- 
tion aperture after the drop has been produced and also 
prevents leaking of the liquid between pulses. 
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FIG. 4: Airflow tube and measurement chamber. 



It is important that the drop generator produce drops 
of constant radius, the primary reason being that the 
size of the drops determines the ratio between Q and 
the measured quantity v x , from Eq. ([!]). A secondary 
reason that we did not initially appreciate has to do with 
Eqs. H and ||. We set v a i r close to v z ^erm so that v z 
is much smaller than v z ,ttirm- Thus a small change in r 
leads to a relatively large change in v z . This dispersion 
disrupts the consistent spacing between adjacent drops, 
and a decreased separation b etween drops is undersirable 
for reasons discussed in Sec. 



[VA6 



With a clean and newly tuned drop generator we get 
remarkably uniform drop radii, constant to about ±0.2%. 
The drop generator also ejects in a consistent downward 
vertical direction along the centerline of the airflow tube. 
At 1 Hz operation, the drop ejector exhibits slow drifts in 
its characteristics with time scales of the order of a week. 
These drifts appear as changes in drop size and as desta- 
bilization that manifests as the appearance of satellite 
drops or inconsistent drop production. Typically, these 
effects can be compensated for by small changes in the 
drive parameters or adjustments of the air velocity or 
both. By the end of the first data set, set 1, the drop 
ejector had destabilized to the point where it had to be 
removed from the apparatus for cleaning and refilling. 
Similarly during the taking of data set 2, the drop ejec- 
tor and air velocity required periodic small adjustments. 
The end of set 2 was caused by increasing instability in 



drop production which could not be compensated for. 
We do not know the reason for this behavior. 

In our drop generator the silicone oil drops are pro- 
duced with a spread of charges, \Q\ ranging mostly from 
to about 10 e. A fe w percen t of the drops have larger \Q\. 
As described in Sec. [V A 1 we used drops with \Q\ < 9.5e 



to maintain good precision in the charge measurement. 
We do not know what sets the charge distribution for a 
particular drop generator. But we have the general obser- 
vation that silicone oil gives narrow charge distributions, 
whereas water, mineral oil and most other organic liquids 
give broad charge distributions, with \Q\ values as large 
as several 1000 e or even larger. 



D. Optical system 

Referring to Fig. ||, the stroboscopic light source con- 
sists of a rectangular bank of 20 LED's emitting at 
660 nm. The pulse length was about 100 /is. The lens, a 
135 mm focal length, //ll, photographic enlarging lens, 
images the measurement region onto the face of the CCD 
camera with a magnification of 2.1. 

The rectangular active image area of the CCD camera 
|| is 4.8 mm in the horizontal, that is, x, direction, and 
6.4 mm in the vertical, that is, z, direction. Hence the 
viewing area in physical space is 2.29 mm horizontally 
by 3.05 mm vertically. We remind the reader that the 
electric field is horizontal. The active imaging area is 
an array of 240 horizontal picture elements (pixels) and 
736 vertical picture elements (pixels). We chose this ori- 
entation of the array to maximize the vertical distance, 
maximizing the number of images per drop. 

Given the magnification and pixel density of the CCD, 
one would expect from geometric optics that the shadow 
of a 20 /im diameter drop would cover 2 pixels hor- 
izontally and 5 pixels vertically. The actual observed 
shadow typically covered 3 pixels horizontally and 7 pix- 
els vertically, and had an intensity variation that was 
approximately a two dimensional Gaussian. This can be 
quantitatively described as the convolution of the simple 
shadow predicted by geometric optics with a point spread 
function that is a result of the diffractive effects due to 
the finite aperture of the lens. We do not and should not 
observe diffractive effects caused by the small size of the 
drops. 



E. Airflow tube and measurement chamber 

Figure ^ shows a slightly simplified, dimensioned draw- 
ing of the airflow tube and the measurement chamber. A 
rectangular duct contains the upward flowing air. It is 
8.3 mm wide in the direction of the electric field and 
31.8 mm wide in the direction perpendicular to the elec- 
tric field. The field plates that define the measurement 
chamber are 51 mm high and 28.6 mm wide. The inner 
surfaces of the plates are machined flat and are in the 
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FIG. 5: Schematic of electronic system. The LED and high voltage is synchronized to the 10 Hz clock, while the drop generator 
runs asynchronously. 



same plane as the inner surfaces of the walls of the air- 
flow tube. The optic axis of the optical system passes 
through the transparent side walls of the airflow tube. 

The air velocity is sufficiently small, with a Reynolds 
number on the order of R e = 50 so that the flow is lam- 
inar. The 203 mm length of air flow tube between the 
measurement region and the air inlet allows the air to 
settle into its equilibrium flow pattern. At equilibrium, 
the velocity profile of the air is approximately parabolic 
across the narrow direction of the channel (x axis). 
Across the long axis, the flow has a roughly constant 
central region and falls to zero at the boundaries ||. 



F. Electronics 

All the electronics of the apparatus, Fig. ||, are hard- 
wired to give reliable timing, independent of the opera- 
tion of the computer. A 30 Hz handshaking signal from 
the CCD is divided down to provide a 10 Hz clock that 
synchronizes the LED strobe, the electric field switcher 
and the computer image acquisition. The switching of 
the electric field, which is driven by the clock signal di- 
vided by 4, operates at 2.5 Hz. This results in a cycle 
where two images are acquired with the electric field in 
one direction, and then two images with the electric field 



in the other direction. These relationships between the 
signals is depicted in the timing diagram of Fig. ||. The 
drop generator is driven asynchronously at 1 Hz. 



G. Data acquisition and storage 

Data acquisition was performed by a single desktop 
computer running Linux. The computer was equipped 
with two special components: a digital framegrabber that 
allowed the capture of image data from the camera and 
a general purpose input /output interface board with dig- 
ital I/O and A/D conversion capability. The additional 
inputs allowed the computer to monitor the state of the 
experiment as well as a variety of environmental vari- 
ables. 

All software was custom written in C. Hardware de- 
pendent code was encapsulated into drivers at the kernel 
level, which allowed a guarantee of synchronization of the 
software with apparatus by a combination of hardware 
and software buffering of the data. 

The overall strategy was to acquire data from the ap- 
paratus and write them to files in raw form for later pro- 
cessing off line. Recall that each image frame contains 
736 x 240 pixels, each digitized to 8 bit accuracy. There- 
fore acquiring data at 10 Hz produces a data rate from 
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the CCD camera of about 2 MB/s, much too large to be 
stored. Since an image contains just a few drops, most of 
the pixels in an image have just the background signal, 
which allows the information to be stored using sparse 
As mentioned in Sec. 



storage. 



[ID 



the typical image of 
drop extends over an area of 3 pixels by 7 pixels. In each 
frame the regions containing drops were isolated using a 
thresholding operation. The position of each drop was 
then measured using a simple center of mass algorithm, 
and for each drop only a surrounding region containing 
13 horizontal pixels by 21 vertical pixels of the image is 
written to the output file. 



H. Data collection 

The search was carried out in two sets described in 
Table Q. At the beginning of the Set 1 the drop ejector 
was operated at 0.5 Hz, then at 1 Hz for the remainder 
of Set 1 and all of Set 2. 

TABLE I: Data collection 
Data set Weeks Number of drops Total mass (mg) 



13 
17 



3 377 477 
13 430 167 



12.1 
58.0 



III. CALIBRATION, ERRORS, AND 
MEASUREMENT PRECISION 

A. Electric field and drop radius 



Rewrite Eq. (|l|) in the form 

(—) 

x \ 6irr] ) 



E 



(5) 



Consider nonfractional values of Q = 
0, ±1, ±2, . . . . Then, as shown in Sec. |IVB 
sured values of v x sharply peak at 



ne, n = 
the mea- 



= 0,±1,±2,.... 

From a fit to the center of these peaks, Eq. (|J) gives the 
fitted E/r ratio. 

The electric field strength, E, is calculated from the 
measured voltage across the electric field plates, known to 
3%, and the plate separation, 8.25±0.01 mm. The plates 
are parallel to within 0.1 mrad. Inserting the calculated 
value of E into the fitted E/r ratio, we obtain the drop 
radius r. 

We have two additional checks of the drop radius, one 
from the measurement of the error caused by the Brow- 
nian motion, Sec. 



IIIB 



and the other from the measure- 
ment of the net downward velocity of the drops, v z in 
Eq. (§). 

The drop radius depends to a moderate extent upon 
the size and shape of the voltage pulse applied to the drop 



generator and to a slight extent upon the age and history 
of the drop generator. However over periods of hours the 
average drop radius could be taken to be constant, with 
fluctuations of ±0.2% for individual drops. Since the 
data were analyzed in one-hour-long blocks, the average 
E/r ratio for any given block was known to much better 
accuracy than this. 



B. Brownian motion and drop position 
measurement errors 



The precision of the determination of the drop charge 
depends upon the precision of the measurement of v x . 
Consider the sequence of position measurements Xi of 
the trajectory of a drop. For two consecutive frames, j 
and j — 1, the velocity measurement v x j is given by 



At 



(0) 



Here, At is the time between successive frames, 0.1 s in 
our case. Since At is known with very good precision, 
the error in measuring v x comes from the error in deter- 
mining the Xi of the drop centers, and from Brownian 
motion. Take the error in centroiding to be normally 
distributed with a standard deviation of a c . 

During the time At between any two successive mea- 
surements of the Xi positions, Brownian motion adds a 
random contribution with standard deviation given by 



Ob 



1 kTAt 
Zirrjr 



(7) 



Here k is the Boltzmann constant, T is the absolute tem- 
perature, r\ is the viscosity of air and r is the drop radius. 

The trajectory of uncharged drops, which have no con- 
tribution to their trajectories due to the electric field, can 
thus be written as 



Xj = x + a c j + °bs 

i=l 

j 1,2,... , ^images 



(8) 



with xq set by the initial position of the drop, and where 
the cr Cj i (crh,i) are normally distributed with a std. dev. 
of a c (<7b)- The analysis for charged drops is similar if 
the effect of the electric field is first subtracted from the 
observed data points. It then follows that 



v X jAt = a b j + a. 



c,3 



T c,j-1 



and 



(v. 



x,3 



^c+rf, j = k 
r ,..,.■) ±1- = { -*l \j-k\ = l 
0, otherwise. 



(9) 



(10) 
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Therefore the total error on any given velocity measure- 
ment, a v , is given by a\ = 2a^ + of, and the centroid- 
ing error introduces an anticorrelation with magnitude 
— <j\ in two consecutive velocity measurements, due to 
the shared position measurement. We use the concept 
summarized in Eq. ( |Io|) and the observed distributions 
of the v x j (after removal of the contribution due to the 
alternating electric field) to separate o~ c from ai,. 
We find that averaged over this experiment 



<j c = 0.31 //m, ct(, = 0.47 /xm, 



(11) 



in the measurement region. Compared to the size of 
an individual pixel on the CCD, the centroiding error is 
small, approximately 1 /30 of a pixel. The value of Uf, ob- 
tained provides an independent check on the size of the 
drops, and is consistent with the size determined from 
the terminal velocity and the electric field drift velocity. 

Equation |ll| shows that the Brownian motion error, 
<7b, is about the same magnitude as the error involved 
in finding the drop position, er c . Therefore substantially 
reducing a c through the use of smaller pixels will not by 
itself substantially reduce the error on the charge mea- 
surement, since the Brownian motion error can only be 
reduced by increasing N mages . 

The final charge measurement of a drop is made using 
a single, detailed best fit to the entire observed trajectory 
of the drop, and the final error on the charge measure- 
ment o q is a result of propagating the errors a c and <7b 
through this calculation. 



ment region, or a systematic uncertainty in the radius of 
the drop of the order of 0.3%. We believe that the decel- 
eration is due to the evaporation of the drop as it falls. 
The magnitude of this effect is small enough such that 
it can be neglected in the calculation of v x . As a side 
note, any systematic uncertainties in the radius of the 
drops are ab sorbed by the calibration process described 
in Sec. Ill A , and do not affect the final charge measure- 
ment. 

Similarly, other possible sources of measurement er- 
ror such as apparatus vibrations, optical distortions and 
CCD array distortions, and patch nonuniformities on the 
electric field plates, were negligible. 



IV. DATA ANALYSIS AND RESULTS 

A. Drop selection criteria 

In this section we use q — Q/e, a measure of the drop 
charge in units of the electron charge. We required that 
all drops used in the analysis meet the criteria in Ta- 
ble [n| The criteria are designed to maintain a charge 
measurement accuracy of approximately 0.03 e and to 
reject irregular drops caused by inconsistent operation of 
the drop generator. 



Other sources of errors 



q < 9.5 criterion 



We looked for other sources of errors, but all are negli- 
gible compared to those in Eq. ([n]) . When we developed 
the upward air flow method we thought about the pos- 
sibility that there might be some horizontal air velocity, 
v x ,air in the measurement chamber, contributing an error 
to a v of order v x ^ r x At. By studying a large amount of 
data we found that the distribution of v x .air x At had an 
rms value of 100 nm, and was a fixed property of the mea- 
surement region. For comparison v e x At was of the order 
of 8 /Ltm. These irregularities in v x ^ a i r are probably due to 
residual surface imperfections in the electric field plates. 
Since the irregularities are constant over long periods of 
time, they can be accurately measured and corrected for. 
For this analysis, that was not necessary. 

Another possible source of error would be a nonunifor- 
mity in the electric field in the measurement region giv- 
ing a horizontal gradient, dE/dx. This would produce a 
horizontal force on the drop's induced electric dipole mo- 
ment. This dipole force acts in addition to the QE force. 
We found such a dipole force to be negligible compared 
to the QE force. 

A small, vertical deceleration of the drops as they fall 
through the measurement chamber was observed. This 
amounted to a change of 30 fim./s in the apparent termi- 
nal velocity of the drops as they fell through the measure- 



For any given drop there is an uncertainty in the radius 
of approximately 0.2% which contributes to the relative 
error on q. The absolute error on q thus increases linearly 
with q. Since the absolute error on q must be kept to the 
order of 0.03 e, restricting the data sample to drops with 
q < 9.5 keeps this contribution to less than 0.02 e. The 
overall charge distribution is such that only a few percent 
of the drops have q values outside this range. 



Primarily, this criterion is a measure of Ni, 



of 



the drop. Brownian motion and centroiding acc uracy , 
characterized by a c and <7f, as described in Sec. IIIB, 



limit the accuracy of the charge measurement, a q . For 
any given drop, the number of position measurements, 
Ni ma g es , and the state of the electric field during those 
measurements, in addition to o\, and <x c , determines this 
accuracy. As noted earlier, Ni mages was of the order of 
15. If a drop has an exceptionally large radius or is falling 
too far from the centcrline of the airflow tube, v z will be 



too large and Ni, 



will be too small. 
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3. x 2 criterion 

As mentioned earlier, the final calculation of the charge 
on a drop is done using a fit to the trajectory of the drop. 
It was required that the \ 2 probability of the fit to the 
drop's trajectory be better than 10~ 3 . This rejects a large 
class of rare artifacts based on the statistical likelihood 
that the observed deviations from the fitted trajectory 
could be attributed to the Brownian motion and cen- 
troiding errors. For example, a drop would be rejected if 
it had an anomalous trajectory due to vibrations in the 
apparatus or due to its charge having been changed by 
collision with an ion during measurement. 



4- v z criterion 

The net downward velocity of the drop, v z , depends 
upon the drop radius and the upward air velocity, v a i r , 
Eq. (|^). This criterion insures consistent drop radii 
within the hour long data blocks by requiring 

\v z ,dro P - v z , b iock\ < 0.124 mm/s (12) 

where v z is the measured value for one drop and v Zt u ck 
is the average value of v z for all the drops in the one hour 
data block. Using Eqs. || and ||, taking v air as fixed and 
using an average value for r of 10 /im, this eliminates 
any drops with r different from the nominal value by 
more than about ±0.5%. 

Recall that the air velocity is approximately parabolic 
and that close to the centerline it is given by 

Vair{x) = Vairfi [l ~ {x/x w ) 2 ] (13) 

where x is the distance along the x axis from the center- 
line of the airflow tube, x w = 4.15 mm is the distance to 
the wall of the tube, and u a ir,o is the air velocity along 
the centerline. Therefore this criterion indirectly restricts 
how far the drop can be from the centerline. 



5. x deviation criterion 

This criterion 

\x - xuock\ < 0.19 mm (14) 

provides a direct constraint on how far a drop may devi- 
ate from the centerline in the x direction. Here Xbiock is 
the average value of x for all the drops in the one hour 
data block. The purpose of this criterion is to elimi- 
nate drops that were produced irregularly. The 0.19 mm 
upper limit in Eq. ( |14| ) was determined by examining 
the distribution of x positions of drops produced during 
normal operation of the drop generation and setting the 
upper limit to eliminate the tails. 
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FIG. 6: The viscous coupling between a moving drop 1 on 
a neighboring drop 2 in still air. The small arrows show the 
vector velocity of the disturbed air. Note that there is a slight 
disturbance at the position of drop 2 that will affect the tra- 
jectory of drop 2. 

6. Minimum distance R between any two drops criterion 

The drops interact with one another through their 
induced electric dipole moments and viscous coupling 
through the air. Consider two drops of radius r, drop 1 
moving with a velocity v x due to the force of the electric 
field on its charge, as shown in Fig. ^. This motion will 
move the surrounding air. At the position of drop 2, the 
velocity of the air in the x direction, V x disturbed air' ^ s 
given by 

^disturbed air = |^rA + cos ^)- ( 15 ) 

Since drop 2 sits in this disturbed air, its v x due to the 
force of the electric field on its charge will have super- 
imposed upon it V disturbed air • This will distort the 
charge measurement. Therefore V disturbed air mus t be 
kept small by keeping R, the distance between the drops, 
much larger than r, the radius of the drop. A large sep- 
aration also serves to minimize the interaction between 
the induced electric dipole moments of the drops, which 
increases as the inverse fourth power of the separation. 
We require 

R > 0.62 mm (16) 

separation between any two drops, which limits these 
forces to a small fraction of QE. 

7. Summary and magnitude of drop selection criteria 

Table || gives the percent of drops removed by each 
criterion averaged over each of the two data sets. The 
total percent of drops removed is also given. Since the 
same drop may be removed by several criteria, the total 
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TABLE II: Drop selection criteria. The entries are the per- 
cent removed by each criterion separately. The bottom row 
gives the total percent of drops removed by all criteria. Since 
the same drop may be removed by several criteria, the total 
percent removed is not the sum of the percent removed by 
the individual criteria. 



Criterion 


Set 1 


Set 2 




<7 


0.4 


0.7 




ffg 


6.0 


0.3 


10 J - 


x 2 


4.7 


2.4 




V z 


4.6 


1.2 




X 


9.4 


5.1 


10 2 r 


R 


12.5 


3.8 




Total 


22.3 


8.7 





-4 -2 2 4 



FIG. 7: The q charge distribution in units of e. 



percent removed is not the sum of the percent removed 
by the individual criteria. 




-0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 



FIG. 8: The q c charge distribution in units of e. 




B. Results 



FIG. 9: The q c residual charge distribution in units of e. 



After the application of these criteria we had a fi- 
nal data sample of 1.7 x 10 7 drops of average diameter 
20.6 A™. The total mass of the sample was 70.1 mg. Fig- 
ure [7j shows the charge distribution in units of e. (The 
asymmetry of the charge distributio n is a property of the 
drop generator as discussed in Sec. IIC.) We see sharp 
peaks at integer numbers of charges and no drops fur- 
ther than 0.15 e from the nearest integer. We emphasize 
that there is no background subtraction here, this is all 



the data after the application of the criteria previously 
discussed. 

To show the shape of the peaks at integer values of q we 
superimpose them in Fig. || using the charge distribution, 
q c , defined by q c = q — N c where iV c is the signed integer 
closest to q. The peaks have a Gaussian distribution 
with a standard deviation of 0.021 e. The absence of 
non- Gaussian tails is what allows this search method to 
be so powerful. 
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In Fig. H we superimpose the valleys between the peaks 
using the residual charge distribution, q r , defined by q r — 
q — Ni where Ni is the largest integer less than q. We did 
not find any drops with residual charge between 0.15 e 
and 0.85 e. In this residual charge range there are fewer 
than 1.17 x 10 ~ 22 fractional charge particles per nucleon 
in silicone oil with 95% confidence. 

Thus this 70.1 mg search did not confirm the one un- 
usual aspect of our previous 17 mg search, Halyo et a/.[|), 
where we found 1 drop with a q r of about 0.29 e. No such 
charge was found in this search. While it is of course still 
possible that the fractional charge found in the 17 mg 
experiment was real, we are inclined to believe that the 
17 mg experiment had a very small background that has 
been eliminated by the improved method of this experi- 
ment. 



CONCLUSIONS 



Comparison with other fractional charge 
searches in bulk material 



In Table III we compare this search with previous, 



larger sample, searches for fractional charge particles in 
bulk matter. No evidence for fractional char ge p articles 
was found in the searches by Marinelli et al. Jll | , Smith 
et ai.[jn|, and Jones et aZ. JH|, similar to the null result 
in the present search. 

In their superconducting levitometer search in nio- 
bium, LaRue et aZ.]l0| claimed to have fractional charge 
particles with e/3 and 2e/3. But Smith et al. Jl^] who 
also searched in niobium using a ferromagnetic levitome- 
ter method did not find any such evidence in a four times 
larger sample. At present the results of LaRue et al. |l(]] 
are not understood and are generally not accepted. 

Our search is by far the largest to date and has the 
smallest upper limit of any search on the concentration 
of fractional charge particles in bulk matter. But it is im- 
portant not to generalize our limit to other kinds of bulk 
matter for several reasons. First, we do not know what 
happens to fractional charge particles that are in natural 
matter when that matter is processed. Note that except 
for the search in meteoritic material by Jones et al. [p.3|, 
all the material in Tabic III is processed. 



Second, if we assume the existence of stable, fractional 
charge particles, we do not know what natural materials 
are most likely to have a detectable concentration. Our 
own thoughts are that the most promising natural ma- 
terial is that found in carbonaceous chondrite asteroids, 
since they are representative of the primordial composi- 
tion of the solar system, having not undergone any geo- 
chemical or biochemical processes. Hence similar to the 
motivation of Jones et al. ]l| , our next search will be in 
meteoritic material from an asteroid. 



TABLE III: Searches for fractional charge particles in ordi- 
nary matter. All experimenters reported null results except 
LaRue et ctZ. JlOj] . There are 6.4 x 10 20 nucleons in a milligram. 



Method Experiment 



Material Mass(mg) 



levitometer LaRue et al. 1 10 1 niobium 1.1 

levitometer Marinelli et al. ll) iron 3.7 

levitometer Smith et al. |12j niobium 4.9 

levitometer Jones et al. [[L3J meteorite 2.8 

liquid drop Halyo et al.M silicone oil 17.4 

liquid drop this search silicone oil 70.1 



B. Remarks on further use of this new method 

The purpose of the new method 0] used in this exper- 
iment was to allow large drops to be used compared to 
the classical method, thus increasing the rate at which 
we could search through a sample and also enabling the 
use of suspensions of more interesting materials. We have 
succeeded in doing this, using drops of about 20 fim di- 
ameter compared to the approximately 10 fim diameter 
used in Halyo et fflZ.Q. In the Appendix we discuss fur- 
ther increasing the search rate by using still larger drops 
and by using multiple columns of drops to increase the 
total rate of drop production. We find that with this new 
method the mass per second search rate can be further 
increased by a factor of the order of 10, but probably not 
by a factor of a 100. 



APPENDIX A: 



INCREASING THE SEARCH 
RATE 



Three are several ways in which the mass per second 
search rate can be increased in this experimental method. 



1. Use of larger drops 

The first way to increase the search rate is to use larger 
drops. Maintenance of the precision of the charge mea- 
surement requires that Ni mages increase in proportion to 



the drop radius. An increase in Ni mages can be accom- 
plished by some combination of a decrease in v z and an 
increase in the vertical length Z, Eq. (Q). However a 
significant decrease in v z requires too fine a balance be- 
tween v a i r and r 2 , Eq. (^|). If we keep v z constant, an 
increase of Z can be attained by an increase in the num- 
ber of vertical direction pixels in the CCD array of the 
camera. Existing CCD cameras with 10 Hz frame read- 
ing rates have twice the 736 vertical pixels used in the 
present camera and larger arrays will probably be avail- 
able in the future. Therefore based on this consideration 
alone, drop diameters of several times 20 /im are feasible. 
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However, there are two problems that must be consid- 
ered for drop diameters larger than 30 to 40 fim. The 
dipole force on a drop in a non-uniform electric held is 
proportional to the third power of the drop dia meter . 
This force was negligible in this experiment, Sec. Ill C , 
but would have to be considered for much larger drops. 
The other problem is that the maintenance of a small 
and constant v z , Eq. (^), requires v a i r to increase as the 
square of the drop diameter, possibly leading to non- 
laminar flow. Therefore without more design and exper- 
imental studies, our conservative conclusion is that the 
drop diameter is limited to about 30 /xm. This would 
lead to an increase of the mass search rate by a factor of 
3 compared to the 20 /mi drops used in this experiment. 



!. Increase of drop production rate per column of 
drops 

Let the drop production rate for a column of drops be 
n per second. Then the vertical separation betwee n drops 
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in a column is R = v air /n. The criteria in Sec 
require R > 0.62 mm. Using v a i r = 2.0 mm/s, this gives 
an upper limit on n of about 3 H z. However our experi- 
ence in this experiment, Sec. II C , strongly suggests that 
a maximum 1 Hz rate is conservative practice, because 
of irregularities in drop production. 



3. Increase in the number of drop columns 



horizontal direction in the CCD array must be increased 
and the space between the electric field plates must be 
increased. The latter requirement means the alternating 
potential applied across the plates must also be increased 
to keep the electric field constant. 

Existing 10 Hz frame rate CCD cameras limit the num- 
ber of columns to three but improvements in these cam- 
eras would probably allow five columns. The correspond- 
ing increase in the electric plate spacing and the potential 
difference is straightforward. 



4. Correction for drop to drop interactions 

It is clear that the primary constraint limiting the den- 
sity of drops achievable in the measurement chamber is 
that of Eq. ([Lsl). To reiterate, interaction between the 
drops due to their induced electric dipole moment and 
viscous coupling requires that there be a minimum sep- 
aration allowable between drops. In the limit that these 
interactions are small, both of these effects can be cal- 
culated from first principles, for example as in Eq. Jig). 
In principle then, it should be possible to subtract the 
effect of these perturbing forces from the measured tra- 
jectory of each drop. Given this, it would be possible to 
relax the constraint on R. Since this possibility requires 
further study, it is not clear to what extent R can be 
reduced and throughput increased. 



In this experiment we used one column of drops, how- 
ever the extension to many columns of drops is straight- 
forward. Of course the horizontal separation between ad- 
jacent columns must meet the requirements of Eq. (|l6|), a 
nominal separation of 1 mm is useful for design purposes. 
The use of multiple columns requires two changes in the 
experimental design, namely the number of pixels in the 



5. Summary 

Putting these estimates together we can see how to 
achieve an improvement on the order of 10 times the 
present mass per second search rate using existing CCD 
cameras. Future cameras will probably allow a factor of 
15 improvement. 
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